A heavier gluino from t-b-r Yukawa- unified SUSY 
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^ Abstract 

Supersymmetric models with t—b—r Yukawa coupling unification and unified gaugino 
masses at the GUT scale- with fi > 0- show a mild preference for light gluino masses 
nig < 500 GeV. This range of rrig is now essentially ruled out by LHC searches. We 
show that a heavier gluino with m- g ~ 0.5—3 TeV can also be compatible with excellent 
t — b — r Yukawa coupling unification in supersymmetric models with non-universal 
Higgs masses (NUHM2). The gluino in such models is the lightest colored sparticle, 
while the squark sector displays an inverted mass hierarchy with rrig ~ 5 — 20 TeV. 
We present some LHC testable benchmark points for which the lightest Higgs boson 
mass irth — 125 GeV. We also discuss LHC signatures of Yukawa- unified models 
with heavier gluinos. We expect gluino pair production followed by decay to final 



states containing four 6-jets plus four H^-bosons plus missing E T to occur at possibly 
observable rates at LHC. 
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1 Introduction 



Unification at M GUT (~2x 10 16 GeV) of t — b — r Yukawa couplings [1, 2] is largely- 
inspired by the simplest supersymmetric (SUSY) 50(10) or SU(A) C x SU(2)l x 
SU{2) R models [3]. It has become clear in recent years[4, 5, 6, 7, 8, 9, 10, 11] that 
imposing t — b — r Yukawa coupling unification has important consequences for the 
sparticle and higgs mass spectrum of the minimal supersymmetric standard model 
(MSSM). The successful launch of the Large Hadron Collider (LHC) has provided 
important new impetus for these studies[12, 13, 14]. For analogous discussion of b — r 
unification, see Ref. [15]. 

The parameter space of 5*0(10) SUSY GUT models for this investigation is given 

by 

mie, m 2 Hu , m 2 Hd , mi/ 2 , A , tan/3, sign(fx), (1) 

where mis is the unified matter scalar mass, rn 2 Hu and m 2 Hd are the GUT scale Higgs 
soft masses, is the unified gaugino mass, A is the coefficient of the soft super- 
symmetry breaking (SSB) trilinear term, and tan /3 is the ratio of Higgs field vevs. In 
order to allow for an appropriate radiative breaking of electroweak symmetry, the two 
GUT scale Higgs doublet masses must be split [16] according to itl 2 Hu < m 2 H . This 
splitting might arise due to 5*0(10) .D-terms (D-term splitting) or via GUT-scale 
threshold corrections [7] (the Higgs splitting, or HS model). With the MSSM super- 
potential parameter fi > 0, this scenario predicts an inverted scalar mass hierarchy 
(IMH)[17] in the squark sector, wherein third generation squarks have masses in the 
few TeV range, while the first two generations of squarks have mass in the 5-30 TeV 
range [6]. The IMH allows one to reconcile a decoupling solution to the SUSY flavor 
and CP problems with relatively low fine-tuning in the EWSB sector. A successful 
implementation of the IMH scheme requires GUT-scale SSB terms to be related as 
Aq ~ 2m^ ~ 4m^ 6 , with unified third generation Yukawa couplings and the unified 
gaugino mass mi/2 on the low side: typically sub- TeV. 

One particularly important prediction concerns the gluino which turns out to be 
the lightest colored SUSY particle. 1 Since m\i% is favored to be small in comparison 
to m 16 , there is some tendency in t — b — r unified models for rrig < 500 GeV, 
which should be within range of SUSY searches at LHC operating with y/s = 7 
TeV (LHC7)[20]. In this case, due to the large 6-quark Yukawa coupling and the 
inverted squark mass spectrum, gluinos are expected to dominantly decay by 3-body 

1 These predictions are obtained under the assumption that the lightest neutralino is the lightest 
MSSM particle (LMP). This class of Yukawa-unified models tends to predict a thermal neutralino 
relic abundance flLMph 2 >> 1 so that the neutralino LMP in this scenario is not a viable dark 
matter candidate. To overcome this drawback, one proposal[18, 19] is to invoke axion physics and 
arrange for the lightest neutralino to decay before nucleosynthesis into an axino, which now plays 
the role of lightest SUSY particle (LSP). A combination of axions and axinos could then make up 
the dark matter content of the universe. 
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modes such as g — > bbx®, leading to final states at LHC consisting of multiple 6-jets 
+MET[12, 13]. 

In fact, recent searches by the ATLAS experiment with less than 1 fb _1 of data 
already exclude m g < 500 GeV by searching for multijet plus missing Et (MET) plus 
one or more tagged 6-jets[14]. Also, direct searches for gluinos and squarks under 
the assumption of unified gaugino masses by Atlas and CMS (again with ~ 1 fb -1 of 
data) typically exclude m g < 550 — 750 GeV (depending on search techniques) [21, 22]. 
Based on this critical input from experiment, the question arises: Are t — b — r unified 
models now excluded by LHC searches, or can Yukawa-unified solution be found with 
heavier gluinos with mass beyond current LHC reach? And if such solutions are found, 
what is the nature of the SUSY signal which is expected in near future runs of LHC7? 

Our main goal in this paper is to determine whether Yukawa- unified solutions 
with a heavier gluino can exist. In fact, we find numerous solutions, some of which 
are presented as a new set of benchmark points with rrig in the range 0.5 — 3 TeV. In 
these solutions, the gluino retains its position as the lightest colored SUSY particle, 
while squarks remain in the multi-TeV range. Thus, we expect in this class of models 
that LHC searches should focus on gluino pair production. However, for these heavier 
gluino solutions, the g is expected to decay via g — > tbxf or tiy®- After t — > bW and 
Xi ~^ X?^ decays, we expect gluino pair production final states to contain typically 
four 6-jets, four W bosons plus MET. These are in rather sharp contrast with models 
containing rrig < 500 GeV, where multi-6-jets+Mi?T final states are expected, but 
without the numerous on-shell W bosons. 

Note that due to potential threshold corrections which could arise from a vari- 
ety of sources including a more complicated Higgs sector, higher order interaction 
terms, etc., we do not insist on exact (or perfect) unification of the three Yukawa 
couplings. Instead, in this paper Yukawa unification realized at the 10% level (or 
better) is considered to yield an acceptable scenario. In practice, we find solutions 
with heavy (~ 2 — 3 TeV) gluino masses that are associated with Yukawa unification 
at a few percent level. Somewhat lighter gluino masses ( 1-1.5 TeV) are accompanied 
by essentially perfect Yukawa unification! As expected, the squark masses display an 
inverted mass hierarchy, with the lightest (third family) squark masses ranging be- 
tween 1 to 10 TeV. The first two family squarks turn out to be considerably heavier, 
of order 8-28 TeV. In the benchmark points that we highlight in this paper, the mass 
of the SM-like Higgs boson is of order 124-126 GeV, a value which is consistent with 
results from recent ATLAS and CMS Higgs searches[23]. 
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2 Phenomenological constraints and scanning pro- 
cedure 



We employ the ISAJET 7.80 [24] package Isasugra[25] to perform random scans over 
the fundamental parameter space. In this package, the weak scale values of gauge and 
third generation Yukawa couplings are evolved to Mq via the MSSM renormalization 
group equations (RGEs) in the DR regularization scheme. We do not strictly enforce 
the unification condition gs — gi — g<i at Mq, since a few percent deviation from 
unification can be assigned to unknown GUT-scale threshold corrections [26]. The 
deviation between gi = g 2 and g 3 at Mq is no worse than 3 — 4%. For simplicity we 
do not include the Dirac neutrino Yukawa coupling in the RGEs, whose contribution 
is usually small [2 7]. 

The various HS model boundary conditions are imposed at Mq and all the SSB 
parameters, along with the gauge and Yukawa couplings, are evolved back to the weak 
scale M z . In the evaluation of Yukawa couplings, the SUSY threshold corrections [28] 
are taken into account at the common scale Msusy — ^ m t L m t R - The entire parameter 
set is iteratively run between Mz and Mq using the full 2-loop RGEs until a stable 
solution is obtained. To better account for leading-log corrections, one-loop step-beta 
functions are adopted for gauge and Yukawa couplings, and the SSB parameters rrii 
are extracted from RGEs at multiple scales m 8 = m^mi). The RGE-improved 1-loop 
effective potential is minimized at Msusy , which effectively accounts for the leading 
2-loop corrections. Full 1-loop radiative corrections are incorporated for all sparticle 
masses. 

The requirement of radiative electroweak symmetry breaking (REWSB) imposes 
an important theoretical constraint on the parameter space. In order to reconcile 
REWSB with Yukawa unification, the MSSM Higgs soft supersymmetry breaking 
(SSB) masses should be split in such way that m 2 Hd /m 2 Hu > 1.2 at Mq [29]. As 
mentioned above, the MSSM doublets reside in the 10 dimensional representation 
of 5*0(10) GUT for Yukawa unification condition to hold. In the gravity mediated 
supersymmetry breaking scenario [30], the required splitting in the Higgs sector can 
be generated by involving additional Higgs fields [10], or via D-term contributions [31]. 
In our Yukawa-unified SUSY spectrum calculations, the lightest neutralino is always 
turns out to be the LMP. 

We have performed Markov-chain Monte Carlo (MCMC) scans for the following 
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parameter range: 



< m 16 < 30 TeV 
< m Hu < 35 TeV 
< m Hd < 35 TeV 

< m 1/2 < 5 TeV 
30 < tan/3 < 60 

-3 < A /m < 3 (2) 

with jj, > and m t = 173.3 GeV [32]. Note that our results are not too sensitive to 
one or two sigma variation in the value of m t [10]. We use m® R (mz) = 2.83 GeV 
which is hard-coded into IS A JET. 

In scanning the parameter space, we employ the Metropolis-Hastings algorithm 
as described in [33]. The data points collected all satisfy the requirement of REWSB, 
with the neutralino in each case being the LMP. After collecting the data, we impose 
the mass bounds on all the particles[34] and use the IsaTools package [35, 36] and Ref. 
[37] to implement the various phenomenological constraints. We successively apply 
the following experimental constraints on the data that we acquire from Isasugra: 

m h (lightest Higgs mass) > 114.4 GeV [38] 

BR(B 3 n+u~) < 1.1 x 10~ 8 [39] 

2.85 x 10~ 4 < BR(b si) < 4.24 x 10~ 4 (2a) [40] 

0.15 < ^if^ Tt/r) ^ SSM < 2.41 (3a) 41 

— BR(B u ^tu t ) S M — v ' 11 

As far as the muon anomalous magnetic moment a M is concerned, we require that 
the benchmark points are at least as consistent with the data as the Standard Model 
is. For a presentation of (g — 2) /t values in NUHM2 models, see [42]. 



3 A heavier gluino from Yukawa- unified SUSY 

In order to quantify Yukawa coupling unification, we define the quantity R t b T as 

_max(y t ,y b ,y T ) 
mm(y t , y b , y T ) 

In Fig. 1 we plot R thr versus the various 5*0(10) model input parameters. Gray points 
are consistent with REWSB and neutralino LSP. Orange points satisfy the mass 
bounds (including mh in the range 115 — 131 GeV and and nig > 0.5 TeV), constraints 
from BR,(B S — > fi + fi~), BR(B U — > rv T ) and BR(b — > S7). Blue point solutions belong 
to a subset of orange points and represent rrih in the range 123 — 127 GeV. In Fig. 
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Figure 1: Plots in m 16 - R tbr , mi/ 2 - R t b T , tan /3 - R tbr , A /m 16 - R tbr , and m~ g - R thr 
planes. Gray points are consistent with REWSB and neutralino LSP. Orange points 
satisfy mass bounds (including rrih in the range 115 — 131 GeV and m- g > 0.5 TeV), 
constraints from BR(B S — > BR(B U — > tu t ) and BR{b — > 57). Blue point 

solutions belong to a subset of orange points and represent m h in the range 123 — 
127 GeV. 
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la), we see, as is well known, that m 16 > 10 TeV for solutions with R tbr < 1.1, as 
required by the inverted scalar mass hierarchy. In Fig. lb), we also see that low 
values of mi/2 are favored. While previous works favored mi/2 < 0.1 — 0.2 TeV, here 
our dedicated MCMC scans show t — b — r unified solutions can also occur for m 1 / 2 
values in the 0.3 — 1 TeV range. Fig. lc) shows that tan/3 ~ 50 — 60 is required, 
while Fig. Id) shows that A ~ — 2m 16 is required for the IMH. Our key result here 
occurs in Fig. le), where we plot the value of m g vs. Rtbr- Here, we find that while 
many solutions occur with m g < 0.5 TeV, there also exist solutions with near perfect 
Yukawa unification with substantially heavier gluino masses ranging up to m g ~ 1.4 
TeV. And if we only require Rtb T ^ 1-1, then some solutions can occur with m g as 
large as 3 TeV! 

4 Heavier gluino benchmark points and implica- 
tions for SUSY searches at LHC 

In Table 1, we list four benchmark (BM) Yukawa-unified solutions from Isajet 7.80 
with trig > 500 GeV. Each BM point also has irih = 125 ± 2 GeV, so all are consistent 
with the Atlas/CMS hint of a Higgs signal around 125 GeV. 

For point 1, with itliq ~ 21 TeV, all the squarks and sleptons are far beyond 
the reach of LHC. However, for this point, m g = 750 GeV, and so gluinos would be 
pair-produced at LHC7 with a cross section of ~ 60 fb[12] . In Ref. [12], LHC search 
strategies assumed a much lighter gluino of mass ~ 0.3 — 0.6 TeV, in which case 
gluino three body decays to bbxl are dominant, and the search strategy was to look 
for collider events containing multiple 6-jets +MET. For point 1, at the bottom of 
the Table we list the dominant gluino branching fractions. In this case, with m g ~ 750 
GeV, the decay modes g —> tix® and g —> tbxj occur at substantial rates: in this 
case ~ 60%. Here, Xi ~^ X^\W a ^ 100% branching fraction, while t — > bW also at 
100%. Thus, gluino pair production for Yukawa-unified benchmarks and a heavier 
gluino lead to final states including four 6-jets, four on-shell W^-bosons +MET. Since 
the Ws decay into hard isolated leptons over 20% of the time, these gluino pair 
production events will contain high multiplicities of isolated leptons, including same 
sign (SS) and opposite-sign (OS) pairs, trileptons and four-leptons! There are few 
SM background (BG) processes that can lead to events containing for instance four 
6-jets plus four isolated leptons. The major BG process would likely be four top 
production: pp — > titiX. 

In addition, xt X2 production can occur at large rates for the Yukawa- unified 
BM points[12]. For all cases listed, xt ~^ XiW a t ~ 100% branching fraction, and 
X2 xV 1 with typically a branching fraction > 90%. Recently, it has been pointed 
out [43] that the process pp — > X1X2 ~^ WhxiXi should be visible at LHC with 



7 





Point 1 


Point 2 


Point 3 


Point 4 




21370 


20230 


18640 


26130 


mi/2 


93.41 


364 


579 


1021 


A /m 16 


-2.43 


-2.13 


-2.09 


-2.11 


tan f3 


57.2 


51 


50 


52 


m H d 


22500.0 


26770 


24430 


34210 


m H u 


13310.0 


23260 


21780 


30590 


m h 


126.7 


125 


124 


124 


m H 


noon 
9389 


o i no 
3192 


1 A K 

3145 


4U66 


m,A 


Oqno 

yozo 


1 1 71 


Q1 OK 

olzo 


/in/in 
4U4U 


m H ± 


yoyu 


^1 Q'? 


11 A 7 
Ol4 I 


4U0 ( 


m- g 


i OU 


1 Q7£ 


loOo 


9nm 
zyyi 


771 -o 

Xl,2 


122, 285 


232, 491 


323,661 


557,1114 


771^0 

A3, 4 


lyzyo, ryzyo 


fiD/18 fin/18 
0U4o,0U4o 


A K7C\ A c;71 
40 I U,40 I 1 


Ooio,ooio 


m Y ± 


zoo, lyoou 


4yo,ouzi 


AKAO 
004,404Z 


1118 fi97^ 




21389,21132 


20230,20115 


18653,18574 


26187,26079 




7389,8175 


3465,5356 


3089,5447 


4376,7901 


m d L ,n 


21389,21513 


20230,20333 


18653,18742 


26187,26304 




7836,8234 


5417,6047 


5534,6584 


8038,9652 


m 0l 


21196 


20128 


18565 


26037 


77lo 3 


15502 


15066 


14032 


19441 


m e L ,R 


21193,21717 


20123,20416 


18559,18779 


26027,26319 


77lf 12 


7490,15463 


8048,15079 


7796,14042 


9984,19455 


QcDAlh 2 


12642 


190 


972 


1377 


Rtbr 


1.06 


1.00 


1.05 


1.07 


BF(~g -> bbxl) 


0.33 


0.13 


0.07 


0.06 


BF(g -> 


0.15 


0.15 
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Table 1: Sparticle and Higgs masses (in GeV). All of these benchmark points sat- 
isfy the various constraints mentioned in Section 2 and are compatible with Yukawa 
unification. Point 1 exhibits a solution near the current reach limit of LHC. Point 
2 exhibits 'perfect' Yukawa unification. Point 3 displays an example of a relatively 
heavy gluino within reach of LHC 14. Point 4 represents a solution with the heav- 
iest gluino (~ 3 TeV) we have in our scans; it is likely beyond reach of LHC. The 
uncertainty in the Higgs mass (m/J estimates is about ±2 GeV. 

a/s = 14 TeV and ~ 100 — 1000 fb _1 of integrated luminosity. This latter gaugino pair 
production signal offers a second corroborative channel for claiming a SUSY discovery 
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in models with lighter gauginos and decoupled squarks and sleptons. In addition, 
in the Wh channel, the Pr{h) distribution may allow a chargino/neutralino mass 
extraction provided a very large data sample is acquired. Likewise, for nig ~ 500 — 800 
GeV, then x° — ► Xi^. In this case, X1X2 production will yield WZ + MET events, 
for which the WZ — » 3£ and possibly WZ — > £ + £~ + jets signatures may be visible 
at LHC7[44]. 

For point 2 in Table 1, we show a case with essentially perfect Yukawa unification, 
Rtbr = 1-0, but with rrtg = 1375 GeV. In this case, the combined branching fraction 
for gluinos into top quark final states has increased to ~ 86%. With such a heavy 
gluino, this case would likely be beyond the reach of LHC7[20]. However, it should 
be within reach of LHC with y/s = 14 TeV (LHC14), which should start operating 
around 2015. Benchmark point 3 in Table 1 shows a case with rrtg = 1853 GeV 
and decoupled scalars. This case, with such a heavy gluino mass, lies right around 
the ultimate reach of LHC14 with 100 fb _1 , in a search for gluino pair production. 
However, the Wh search channel from X1X2 production may be competitive with 
gluino pair searches in this case assuming 100-1000 fb _1 of integrated luminosity. 

The last point 4 in Table 1 corresponds to R t b T = 1-07, but with nig = 2991 GeV 
and decoupled scalars. This is the case with the largest gluino mass we were able to 
find while requiring Rtbr < 1.1 and rrih ~ 125 GeV. This case would likely lie beyond 
reach of LHC14 for any luminosity upgrade. Detection of a SUSY signal in this case 
would likely require a pp collider with -y/s ~ 40 — 100 TeV. 

5 Conclusion 

Previous papers examining t — b — r Yukawa-unified models with gaugino mass unifi- 
cation and /i > have focused on solutions with rather light gluinos: nig < 0.5 TeV. 
These models are now likely all excluded by recent or soon-to-be-released LHC SUSY 
searches. In light of these earlier results, we were motivated to examine if Yukawa- 
unified solutions with heavier gluinos could exist, while also requiring nih ~ 125 
GeV, as is recently hinted at by Atlas and CMS. Using dedicated MCMC scans over 
S'O(IO) parameter space, we have found solutions with excellent Yukawa unification 
and rrtg ranging up to 1.4 TeV, well beyond current LHC search limits. Loosening 
the Yukawa- unification criteria to Rtb T < 1-1, we even find solutions with rrig nearly 
3 TeV. 

We have listed four 5*0(10) benchmark points with m g spanning the range 0.75 — 
2.9 TeV. Regarding LHC SUSY searches, we note that these heavier gluino solu- 
tions will be characterized by gluino pair production at LHC, followed by decays to 
final states including four fe-jets, four on-shell W bosons +MET. The gluino pair 
events should be rich in multiple isolated leptons plus 6-jets, and the dominant SM 
background will likely arise from four top production. 
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